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The space-like electromagnetic form factor is expressed in terms of the overlap integral 
of the initial and final meson wave functions written as Lorentz covariant distributions 
of internal momenta. The meson constituents are assumed to be valence quarks and an 
effective vacuum-like field. The momentum of the latter represents a relativistic general- 
ization of the potential energy of the quark system. 

The calculation is fully Lorentz covariant and the form factors of the charged mesons 
are normalized to unity at t—0. The numerical results have been obtained by freezing 
the transverse degrees of freedom. We found r^=0.434 fm^, =0.333 fm'^, r^o =-0.069 
fm^ by taking m„,d=430 MeV, ms=700 MeV. 
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I. INTRODUCTION 

In the standard model and the lowest order of perturbation the elastic electron meson scattering is the 
result of one photon exchange between the electron and one of the quarks in the meson. Formally the 
hadronic matrix element entering the expression of the scattering amplitude is parametrized in a model 
independent Lorentz covariant way as: 

< Af(P')l t/.(+(X3,0) J^^iO) C/s(0,-oo) \M{P) > = = Unit) [P + P'r (1) 

where is the elementary electromagnetic current expressed in terms of free quark fields, ?7s(t, t') is 
the time evolution operator of the system due to quantum fluctuations of the colour field and t = [P—P'Y 
is a space-like momentum transfer. All the information concerning the quark structure of the meson and 
the effects of the quantum fluctuations of the coloured field is contained in the electromagnetic form factor 
fem(t). This makes the calculation of the electromagnetic form factor a rather complicate problem. 

As shown in a careful analysis [Q, perturbative QCD calculations solely cannot explain the observed 
pion and nucleon form factors at presently accessible values of i . The dominant contribution seems to 
come from the hadronic wave functions which are essentially nonperturbative. 

In fact this was to be expected from the uncertainty principle of quantum mechanics. Point-like quarks 
and a real local interaction between them could be observed only asymptotically. At low momentum 
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transfer only dressed quarks, the so-called constituent quarks, are observable. Moreover, the electromag- 
netic current expressed in terms of constituent quark fields is only approximately local and the binding 
may also induce sizable effects. 

In the present paper we calculate the form factors at low momentum transfer, where mesons exhibit 
a stable, time invariant structure in terms of constituent quarks. We assume that, to a good extent, 
the electromagnetic current expressed as a product of constituent quark fields may be considered local. 
Consequently the matrix element (|l|) will be related to the overlap integral between the wave function 
of the final meson and that of the initial meson where one of the quarks has been replaced by the quark 
interacting with the electromagnetic field. We therefore consider that the space-like form factor indicates 
to what extent the initial system, where one of the quarks has been replaced by the recoiling quark, 
represents the final meson. Similar significance can be assigned also to the semileptonic form factors. 

In this framework, at low space-like t the form factor fem {t) is the effect of a stable meson structure [|l| 
and one would expect it to decrease continuously with \t\ because the larger is the momentum transfer, 
the less probable is to find the recoiling quark in the final meson. The decrease is then nothing else but 
the consequence of the worse and worse matching of the initial with final wave functions. 

Then the calculation of the form factors at low momentum transfer amounts in finding a Lorentz 
covariant internal wave function for the meson as a bound system made out of independent constituents. 

One of the best known solution to this problem is the Bethe Salpeter wave function. Completed with 
the expression of the dressed quark-photon vertex satisfying the Dyson-Schwinger equation it led to 
remarkable results for low energy observables |^ , [Q . It must be noticed however that the formalism has 
ambiguities related to the definition of the relative quark momentum and that the presence of negative 
energy states in the quark propagators violates the real quark content of mesons. Also, in order to satisfy 
the usual normalization condition of the form factor additional assumptions seem to be necessary [0,0. 

Quark models also used solutions of the Schrodinger equation with suitably chosen potential terms. 
The calculations have been performed at i = and the decrease with |t| has been introduced from 
phenomenological reasons 

Results independent of the specific form of the confining potential have been obtained in the light-front 
approach where the potential is formally included into the free mass term of the quark-antiquark 
system. The form factors can be calculated at any i, but manifest Lorentz covariance is lost because of 
the difference between "good" and "bad" components of the current. 

In this paper we present an original way to take into account the confinement Q. The specific as- 
sumption of the model is that at low energy the hadrons look like stable systems made of valence quarks 
and of an effective vacuum-like field $ representing countless excitations of the quark-giuonic field which 
cannot be observed one by one. The sum of all their momenta is the 4-momentum of the field which 
is supposed to be the relativistic generalization of the potential energy of the quark system. 

The main reason for introducing this effective description is suggested by the exemples taken from 
relativistic field theory, where a bound state is the result of an infinite series of elementary interactions, 
but not the instantaneous effect of an elementary process |^ . 

We therefore conjecture that the binding forces are the result of a time average of elementary processes 
which cannot be seen one by one if the time of observation is rather long (i.e. if the momentum transfer 
from the projectile to the target is rather low). We consider that the countless elementary excitations 
continuously taking place in a bound system deserve an average description in terms of an effective field 

not an individual, perturbative treatment. Furthermore, the presence besides the valence quarks of an 
effective field $ whose 4-momentum does not satisfy any mass constraint may also be useful for ensuring 
a definite mass to the compound system. Indeed, as it is well known, a system made of on-mass-shell 
particles having a continuous distribution of relative momenta does not behave like a single particle 
because it does not have a definite mass |^. An effective field may improve the situation by playing the 
role of a potential whose contribution adds to the quark kinetic energy yielding a definite value for the 
meson mass. 

The main features of the model are presented in the next section. We give the generic form of the meson 
state and calculate the expression of the electromagnetic form factor by using the canonical formalism of 
field theory. 
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The third section contains the numerical results obtained with frozen transverse degrees of freedom. 
In the last section we draw some conclusions concerning the applicability of the model and comment 
on the values of the parameters yielding the best fit. 



II. THE MODEL 



The main feature of the model is to suggest the generic wave function of the bound quark system 
representing the meson. A first requirement for the wave function is to be a Lorentz covariant function 
of the independent quark coordinates because the quarks really behave like independent particles in the 
interaction with the electromagnetic field. It is obvious that in order to satisfy this requirement without 
introducing unphysical coordinates it is necessary to work in momentum space. This will also help 
expressing easily the mass shell constraints and the conservation laws. For this reason the internal wave 
function of the meson will be given under the form of the distribution function of the internal momenta. 

Another constraint for the meson wave function follows from the requirement to satisfy a single particle 
normalization condition. Accordingly, the integral over the internal momenta must converge, the relative 
momentum cut-off being provided by the internal distribution function of the quark momenta. 

Recalling our conjecture that the binding potential can be described with the aid of an effective field 
we assume that the generic form of a single meson state is fl]: 



where a+, 6+ arc the creation operators of the valence qq pair; u, v are Dirac spinors and Tm is a Dirac 
matrix ensuring the relativistic coupling of the quark spins. In the case of pseudoscalar mesons rj\/ — ij^. 
The quark creation and annihilation operators satisfy canonical commutation relations and commute with 
$+((5), which represents the mean result of the elementary quark-gluon excitations responsible for the 
binding. Their total momentum is not subject to any mass shell constraint and, in some sense, it is 
just what one needs to be added to the quark momenta in order to obtain the real meson momentum. 
This is in agreement with our assumption that Q is the relativistic generalization of the potential energy 
but we shall avoid for a while making any definite assumptions about the momentum carried by the 
effective field. 

The internal function of the meson is the Lorentz invariant momentum distribution function ip{p, q; Q) 
which is supposed to be time independent, because it describes an equilibrium situation. This means 
that it does not change under the action of internal strong forces and hence the time evolution operators 
Us{t,t') in eq. ([^) can be replaced by unity. As mentioned above, the main role of is to ensure the 
single particle behaviour of the whole system, by cutting off the large relative momenta. 

In the evaluation of the matrix element (|l|) we shall use the cannonical commutation relations of the 
quark operators 



and the expression of the vacuum expectation value of the effective field which is defined as follows 



(0| $ (Qi) $+(g2) |0) = (yTo)-i Xe^ iQ2-Q.). ^ ^^^^i (^yj^^yi ^(4)(q^ _ 



where V is the meson volume and Tq is the characteristic time involved in the definition of the mean field 
$. To is the time needed by the bound state to be formed and hence we expect it to be of the order of 
the hadronization time. 




xS^^\p + q + Q~P) ip{p,q;Q) X uip)TMviq) '^HQ) a\p)b\q)\{)) 



(2) 



{adk),a]{q)} = {bdk),b]{q)} = (2^)^^ 5.,,5^^Hk - q) 

J J Tin 



(3) 
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It is important to remark that the definition (||) is compatible with the norm of the vacuum state if 
one takes $(0) — 1. We notice also that the relation has the character of a conservation law, just 
like the commutation relations (||) , both of them being necessary for the fullfilement of the overall energy 
momentum conservation in the process. 

As a first test of the model we evaluate the norm of the single meson state (|^) according to the usual 
procedure. The exponent in the integral Jrp cIXq e'^i^iP)-^'^^ ))^o coming from the j''*-' functions in eqs. 

and (||) shall be put equal to and the integral equal to T, because the uncertainty in the meson 
mass is expected to be much smaller than l/T. 

Observing that T is the characteristic time involved in the definition of the effective field $ in the case 
of a moving meson, we write it as T = Tq and get: 

{M[P')\M{P)) =2S (27r)3 5(3)(P-P') J (5) 

where 

J = ^ [d^P^d^^ d^Q 6i^Hp + ^ + Q- pMp, Q)P t. f ^ i-^ = 1. 

2M V J Bp Cq \ 2mi 2m2 / 

(6) 

This a remarkable result because it shows that the wave function of the many particle system repre- 
senting the meson can be normalized like that of a single particle if the integral converges. 

As a matter of consistency, we also remark the disappearance of the rather arbitrary time constant Tq 
from the expression of the norm. 

We evaluate now the matrix element (Q) proceeding in the same manner as before. The expression of 
the electromagnetic current written in terms of free quark fields is 

JenAx) - ^^^Mxh^Mx) (7) 

where is the fraction of the proton charge carried by the quark i. Introducing (|^) between the meson 
states (|^) and using the relations and to eliminate some integrals over the internal momenta we 
obtain after a straightforward calculation: 



where 

=«^j(4)(^ + ^ + Q_p/) ^W(^k,q;Q) Tr [75(fc + mi)7^(p + mi)75(-g + m2)] (9) 
=K2 S^^^p + k + Q - P') <pf{p,k-Q) Tr [75(-^ + m2)7M(-'7 + "^2)75(p + mi)] . (10) 

The two terms in (^) represent the contributions of the valence quarks, k is the momentum of the 
recoiling quark, P' is the final meson momentum and ^pij are the momentum distribution functions of 
the initial and final mesons respectively. 

In the following we shall work in the Breit frame where the momenta of the initial and final mesons 
are P ~ {E, 0, 0, P) and P' = [E, 0, 0, —P) and the electromagnetic form factor reads as: 

= 7mT=t 

In this frame it is an easy matter to show that /em(0)=l in the case of a charged meson. The 
demonstration makes use of 5^'^\k + q + Q) to eliminate the integrals over k in 7^^^^ and of the identity 
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(p + m)7o(p + to) = 2ep [p + m) to reduce the number of projectors. Performing a similar operation on 

the norm, one gets 

To = 2M (ki - K2)J (12) 



(2) 

Tq and proceeding like in the case of the norm, one gets 



which means 

/e,n(0) = 1 (13) 

if the meson wave function is properly normalized. 

In order to calculate the form factor at i ^ we start by using the (^'^'^^ functions to eliminate the 
integrals over the momenta q and k in the expression of T^^^'' and over p and k in the expression of t/i^\ 
After performing the traces over 7 matrices we get 

'^if^'' = / d.(j)p d^Q — S{ep + eq + Qn - E) S{ep - efe) 

X lp^{p,q■,Q) ipf\k,q;Q){qf, t + 2P-Q {kf, - p^) + {k^ + pf,)[£ + ^ 2mi{mi + 1712)]} (14) 



where 



Next, by writing 



£={E-Qof + P^-Q^ + ml-ml. (15) 



^Ji^k - ep) - -^S{p, - P) (16) 



±-Siep ^e,^Qo-E)^^^S (cos ^ - ) (17) 



and 



we perform the integrals over pz and <j)p in eq. ( |14| 
Then Tq^'' becomes: 

-d) - ^1 /■^4n f^w. ..n^ 



4.V1P J J^^^^ 2ptQt V I - COS'' (Pp 

x{2ep [£ - 2TOi(mi + TO2)] + {E - Qo) t}. (18) 

The integration limits over Cp result from the kinematical constraints > to^ + P^ and cos^^p < 1 
which give: 



{E - Qo)£ + Qt^JS^ - mE - QoY - Ql]{ml + P^) 

Gr,M — ~ 

2[(^-Qo)2-gi.] 



_{E- Q„)£ - Qt^£^ - A[{E - Q^f - Ql]{ml + P^) 
""""^ ~ 2[{E-QoY-Ql] ■ ^^^^ 

("2") 

The term 7^ can be proccessed in the same manner, giving a similar expression. 
Using the above results it is possible to calculate the electromagnetic form factors at any any by 
choosing an appropriate function ip. In principle, the calculation does not imply any other approximation, 
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but it is hard to believe that the multiple integral entering the expression of the form factor can be 
performed exactly. 

The expression (|l^) is, of course, valid for t ^ 0, but the infinite value one gets in the limit t seems 
to contradict the normalization of the electric charge ( |l3| ) which has been demonstrated previously. 

This is a disturbing situation which deserves a careful examination. Looking back, we remark that the 
contradiction comes from the evaluation of some S functions: 

S{p + q + Q- P)5{ep + e, + Qo - E{P))5{k + q + Q' - P')6{ek + e, + Qo - E{P'))5^^\Q - Q') (20) 

which have been written as 

5{p + q + Q- P)5{e.p + e, + Qo - E{P))6{p- k - 2P)6{ep - ek)6'-^'> (Q - Q')6{Qo - Q'o) (21) 

at t y^O, while at t—0 they have been written as 

6{p + q + Q)S{ep + eg + Qo- M) x 5{p~ k)5{Qo - Q'o - M + M')d^^\Q -Q')^ f e''(*^-*^''^"VXo. 

27r J 

(22) 

In the last expression the integral has been replaced by T because it was assumed that the uncertainty 
in the meson mass is much smaller than T~^. As it was shown above, at t = the magnitude of T does 
not really matter because it disappears from the expression of the norm. This is not the case at t ^ 
where T remains in the final expression of the form factor. Searching for a way to cure the disagreement 
between the two cases we conjecture that T must be seen as the overlapping time of the initial and final 
systems. Of course, T depends on the reference system we consider. On the other hand, the form factor 
is scalar under Lorentz transformations which means that we have to choose a particular reference frame 
and write T in terms of Lorentz invariant quantities. (We recall that a similar problem occurs in the 
definition of the scattering cross-sections, where the flux of the incident particles is written in terms of 
the relative momentum in the center of mass system.) We choose the Breit frame and consider that T 
is the time necessary for an object of length L to pass along another object of length L. We put then 
T = where L ~ Lq is the Lorentz contracted length of the meson box of size Lq in its rest frame. 

E 

The numerical results quoted in this paper have been obtained by replacing the symmetry scheme based 
on the full Lorentz group with the symmetry under the coUinear group which is equivalent with the flux 
tube model with frozen transverse degrees. In this respect our working approximation may be considered 
as opposite to the light front approach ||] where the transverse degrees of freedom are explicitly taken 
into account, while the confining potential is assumed to give a fixed contribution to the free meson mass. 

In the present approach the relativistic generalization of the potential energy of the quark system is 
explicitly taken into account as the 4-momentum of an effective field and the mean contribution of the 
transverse degrees of freedom to the meson energy is included in the quark masses. The longitudinal and 
the temporal degrees of freedom are then the only active and the multiple integral in eq. ( p^ reduces to 
a simple one. 

In this case the expressions of T^^^'' becomes: 

r^' = U fdQodQ.^p^f^ 5^'Hp + q + Q- P)S^'Hk + q + Q-P') 

VI J 2ep 2eq 2ek 

X(p,(p, q; Q)^^p{k, q; Q) {{cp + eu) [{E - Qof - Ql - P^ - (ml - m2f] + t} (23) 

Proceeding in the same manner as in the three dimensional case we get: 

J_ = —S(p, - P) 
2ek 4P ' 

-^5{ep + eg + Qo-E) = S{ml + Ql-{E-Qo- ep)^) . (24) 
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Introducing now ( p4| ) in ( [2^ ) and performing the replacement E — Qq — ei = m2 a; where ei = y^rn^+P^ 
we obtain finaUy: 



(1) 



1 / TOl \ (1) 



(25) 



In the above relation (pi, tp^^-* are the internal wave functions of the initial and final mesons and the 
integration limits have been deduced from the conditions > or 62 > TO2 which make the square root 
real. 

(2) 

The expression of can be obtained immediately by interchanging toi with m2 and then the space- 
like electromagnetic form factor of pseudoscalar mesons reads: 



fera{i) 



TOl TO2 



dx 



M2 Vo Lo J I V^2^ 



TOl \ (1) / TO2 \ (2) 

Kl ^ —X +11 (fiiipj + K2 I + 1 ) 'Pi'i^f 



(26) 



where ei 2 



TO 



1,2 



p2 . We remark that Vq does not really appear at the end because it will be 



cancelled out by a similar factor coming from the normalization condition of the internal wave functions. 
This is not the case for Lq, which still remains in the expression of the electromagnetic form factor at 
t ^ 0. However, keeping in mind that the form factor is really normalized to unity at t = 0, we shall take 
Lq as a free parameter enabling us to impose this constraint on the expression (|2q). 

A first consequence of the relation (^ ) , independent of the particular form of the internal wave functions 
is that f^„i^'^ {t) = 0, while f^'^ (t) ~ (to^ — to^) and in principle does not vanish. 



III. NUMERICAL RESULTS 



The numerical results we quote below have been obtained for 

fip, q; Q) 



where A/" is a normalization constant. With the notations in eq. (|2^) one has 



(1) 



2 I 



(27) 



(28) 



(TOiTO2)2" 

Observing that for n = 1, 2 the integral in the expression ( |2^ ) can be performed analytically, we find: 



fern (^) 



A'P Vo Lo 

AP Vo Lo 
I 7 

+4 



J2 

i=l,2 

i=l,2 
„2 



2Trmf 



4to^ - i - 4TOf t ^4to2 - t - ^f~t 



(29) 



8mf 



(4to| - t) t2 - 4TOf t 

-4 \ ^4to2 - t + ^J~t 



In- 



^t'^ - im'^t (4to2 - i)^t2 „ 4^2^ ; ^4to2 - t - ^J~t 



(30) 



which gives then: 
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/i^(0) = 



dt 
limt^ 

dt 
linif^ 



M2 Lo 



t=o 



i=l,2 



2 



'*=o AP Vo Lo 



1 



/l^J(t)^_lnH). 



4 
15 



(31) 
(32) 

(33) 
(34) 

(35) 

(36) 



The comparaison of the theoretical results with the experimental data has been done by looking for the 
values of the parameter n and of the quark masses giving the lowest possible x^, where: 



X 



^2 ■ 



(37) 



Here Ui are the experimental errors and the sum is performed over all the experimental points quoted in 
1^ . The values we found for the pion and kaon electromagnetic form factors are given in Figs. 1 (a) and (b) 
respectively. They have been obtained for n=1.5 and the quark masses m„^d=430MeV, ms=700MeV. Our 
fit has xj=69.4 for 45 experimental points and respectively x^=12.2 for 15 experimental points, while the 



simple pole fit has x^=52.5 and x|:= 
the experimental values are {r'^)exp 
fm^ [§. 



:12. We also get r^- 
=0.431±0.01 fm2, ( 



=0.434 fm^, r 



=0.333 fm2, 
=0.34±0.05 and (r^^o) 



2 

K+- 



-0.069 fm^ while 
= -0.054±0.026 
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We notice that the theoretical results are very sensitive with respect to the up and down quark masses. 
This let us believe that the numerical results quoted above may be improved by introducing a difference 
between the lightest quarks. We notice also that the values of the quark masses found from the best fit 
are a bit larger than the currently used constituent masses [0j|^.This is normal taking into account that 
the contribution of the transverse degrees of freedom has been included in the quark masses. 



IV. COMMENTS AND CONCLUSIONS 



The model we presented in this paper is a relativistic model for bound states. It enabled us to calculate 
the electromagnetic form factors by means of the overlap integral over the internal wave functions of the 
initial and final mesons. 

The internal function we used is just a suitable trial function serving to illustrate the qualities of the 
model, not the solution of a dynamical scheme. We remark that in the infinite momentum limit it does 
not have any of the usual asymptotic forms quoted in |^,^. This is because the distribution amplitude 
( p7| ) writes in terms of real Lorentz scalar products which introduce "bad" components of momenta. It 
demonstrates, however, that the existence of the meson form factors can be explained through the stable, 
time invariant structure of the mesons and also, that the colour fluctuations contributing to the quark 
form factor can be neglected at low momentum transfer. 

The remarkable point is that our results arc unexpectedly consistent, as it can be seen from the 
comparison of the pion radius with the box length Lq. First we mention that Lq cannot be infinite or 0. 
In the first case fem{t) = for any t 7^ and mesons would behave like rigid bodies with respect to the 
electromagnetic interaction. In the second case mesons would be structureless. 

In our specific case Lq can be calculated from the normalization conditions of the meson state (||) 
and of the electromagnetic form factor (|l^). Freezing two spatial degrees of freedom, using the internal 
function (p7|) and performing the replacements A/ — Qq = (mi + 7712)0: and Qz = {nii + m2)y we get: 



J = 



2MVo (mi 



m2 



\2n-2 



dx 



dy 



[x^ — 



(mi — m2)^ 
(mi+m2)^ 



1 (38) 



where C = \Jx'^ — \ and the integration limits are the consequence of reality of the norm. 

Then, by using the values of the parameters giving the best fit, we obtain in the pion case Lq ~ 0.46 
fm, not far from the pion radius which is roughly 0.65 fm. 

It is important to notice that the value we get for niay help us to put an upper limit for the time 
To involved in the definition of the effective field <I>. To this end we recall that Lq is related with the 
overlapping time of the systems representing the mesons and conjecture that the model is valid up to 
meson energies where this time is sensibly larger than the time involved in the definition of the effective 
field $. This means ^ p° » ^j^- Then, if we expect the model holds at least up to t = — 1 GeV^ we 
get To « 0.7 IQ-^^s. 

A last point we comment here concerns the continuity of the form factor from the space-like to the 
time-like region. In principle our expression ( p5| ) can be continued at i > 0, but it will have no significance 
in the absence of a dynamical scheme of the strong interaction. 



Acknov^rledgements The author thanks Prof. H. Leutwyler for the kind hospitality at the Institute 
of Theoretical Physics of the University of Bern where this work was completed during author's visit 
supported by the Swiss National Science Foundation under Contract No. 7 IP 051219. The author 
thanks Fl. Stancu, University of Liege, for useful suggestions and careful reading of the manuscript. The 
partial support of the Romanian Academy through the Grant No. 329/1997 is also acknowledged. 



9 



[1] N. Isgur and C. Llellwyn Smith, Phys. Rev. Lett. 52,1080 (1984); ibid. Nucl. Phys. B 317, 526 (1989). 
[2] A. V. Efremov and A. V. Radyushkin, Phys. Lett. B 94, 254 (1980); N. G. Stefanis, W. Schroers, H.-Ch. Kim, 
Phys. Lett. B 449, 299 (1999) ; B. Mehc, B. Ni zic, K. Passek, Phys. Rev. D60, 074004 (1999); V. M. Braun, 



A. Kodjamirian and M. Maul, hep-ph/9907495 



[3] H. Ito, W. W. Buck and F. Gross, Phys. Lett. 52, 28 (1990); ibid. B 287, 23 (1992); ibid. B 351, 24 (1996); 

M. Sawicki and L. Mankiewicz, Phys. Rev D 37, 421 (1988); ibid. 40, 3415 (1989). 
[4] C. D. Roberts, Nud. Phys. A 605, 475 (1996). 

[5] B. Grinstein, M. B. Wise and N. Isgur, Phys. Rev. Lett. 56, 298 (1986); N. Isgur, D. Scora, B. Grinstein and 
M. B. Wise, Phys. Rev. D 39, 799 (1989). 

[6] P. L. Chung, F. Coester and W. N. Polizou, Phys. Lett. B 205, 545 (1988); W. Jaus, Phys. Rev. D 41, 3394 
(1990); ibid., Phys. Rev. D 44, 2851 (1991); F. Schlumpf, Phys. Rev. D 50, 6895 (1994); H.-M. C hoi and 
C.-R. Ji, Nucl. Phys. A 618, 291 (1997); ibid. Phys. Rev. D 59, 034001 (1999); ibtd. \iep-ph 9908431 . 

[7] L. Micu, Phys. Rev. D 55, 4151 (1997). 

[8] S. R. Amendoha et al. (NA7 coIL) Nucl. Phys. B 277, 168 (1986); W. R. Molzen et al. Phys. Rev. Lett. 41, 
1213 (1978). 

[9] Particle Data Group, A. Manohar, Eur. Phys. J. C 3, 337 (1998). 



10 



